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Abstract 

We present the method that the pair production of particles is described as a Fock space state 
vector, where the pair production stems from receiving energy from, as an example, gravitational 
background field. At the same time we show a sort of new mechanism that in the diffeomorphism- 
invariant Dirac sea the Weyl symmetry is broken due to the pair production and this provides 
us a new physical interpretation of the origin of Weyl anomaly. In the present paper at the first 
place, we consider two kinds of background space-time which are connected each other by the Weyl 
transformation. It is shown that the vacuum, i.e. Dirac sea, in the “curved” space-time after 
making Weyl transformation is nothing but the pair produced state provided we consider the state 
in terms of the Fock space basis in the “flat” space-time before making the Weyl transformation. 
When summing up with respect to the time the contribution of pair production of each time-slice 
we obtain the Ricci scalar R and then it is related to the trace part (FJj) of the energy-momentum 
tensor to produce Weyl anomaly. 
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1 Introduction 

It is well known that the virtual pair production and destruction occur always when the external field 
acts to the vacuum. Further these virtual pair produced particles gain the energy from the external 
field and then these particles appear as real particles. In quantum electrodynamics (QED) or similar, 
in particular the case of electron-positron in the presence of the external electric field this is known as 
the Schwinger effect [T] and has been well analyzed theoretically. 

So far the analysis have been done by calculating the vacuum energy density E and then the 
imaginary part ImK is considered as the particle-antiparticle pair production rate defined as 

r = 2Im£'. 

It is physically interpreted that the pair production is due to the tunnel effect, because the energy has 
the imaginary part. 

A direct physical clarification of the pair production in the background field was made by two of 
us of the present authors (H.B.N. and M.N.) [2]. When the electric field is applied to the negative 
energy electrons in the Dirac sea, the right-mover particles and left-mover holes gain the energy and 
are pumped up or down, the pair production occurs. This mechanism is depicted in Figure [1] where 
the blue circle denotes the electron and the white circle is the positron. 

The purpose of the present article is twofold. The first purpose is to reveal in a quantum field 
theoretical method the dynamics of the pair production. We show in more detail and precisely that 
the pair production state is a superposition of all possible particle and antiparticle state. In fact we 
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Figure 1: The pair production due to being pumped up or down from the Dirac sea in [2]. 


will show that the pair production state is indeed a superposition of one pair, two pairs, three pairs, 
... of particle and antiparticle with various energies. 

The second theme of the present paper is a kind of new physical derivation of the Weyl anomaly. 
In 1974, D. M. Capper and M. J. Duff discovered the Weyl anomaly [3] that the trace of the energy- 
momentum tensor for a massless Weyl fermion in a gravitational background is not as classically 
expected to be 0, but rather 


{Tp (X ^R. (1.1) 

Here ^/g is the square root of the determinant of the metric tensor and R is the Ricci scalar. Later 
many recalculations of this anomaly has been made (see the review, for example a). The present 
article is an attempt to obtain a physical understanding of how this anomaly comes about, i.e. to 
reveal the mechanism that the pair production from the Dirac sea leads the Weyl anomaly. As such 
it is also inspired by the interpretation of the chiral anomaly as pumping up or down states from 
or to the Dirac sea [2] of two of us (M.N and H.B.N). An another attempt to understand the 1 + 1 
dimensional Weyl anomaly by our thinking on the Dirac sea we presented in [5]. 

As setting up we consider 1 + 1 dimensional Euclidean cylinder type space-time with the conformal 
flat metric In this space-time we investigate the vacuum-to-vacuum transition 

amplitude of the massless Dirac fermion which has Weyl symmetric action. We then make a basis 
transformation of (or, in other words, pull back) the amplitude into the Fock space in the space-time 
with another metric Then we can show that the Dirac sea in the space-time specified 

by the primed metric g'^^^ ranging from — oo-past to +oo-future is expressed on the basis in the space- 
time with as a pair produced state on each time-slice. By summing up these pair production 
contribution from —oo-past to +oo-future, the difference between the vacuum (or Casimir) energy in 
the space-time of the primed metric g'^^, and that of the unprimed metric g^u gives rise to the Ricci 
scalar. In this way we finally obtain the Weyl anomaly. 

The organization of this paper is as follows. In the following section [2] we summarize the construc¬ 
tion of the Weyl particles and Dirac sea in the gravitational background. In the section [3] we then 
formulate the state description of the pair production and evaluate the Weyl anomaly using the Dirac 
sea formalism. In section 0] we conclude and give an outlook as to how possibly to extend the use of 
our formalism. We have in addition three appendices. 
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2 Complex fermion and the Dirac sea 


In order to describe pair production state as a quantum state in a Fock space, we consider a complex 
fermion 


'F = 


t 


( 2 . 1 ) 


on the 2-dimensional Euclidean cylinder with 

X = ,x^) = (t, x), with — oo < t < -|-oo and 0 < x < 2 ti. 


We adopt the notations that the Roman indices i, j, • • • and Greek indices • • • represent flat and 
curved space-time respectively. The diffeomorphism-invariant action for T reads 

S['^,g] = ^j <fxy/^-^\x)^^^^et{x)V^^{x) 

= ^ / (fx^/g • 7Vei^(^)^(®)} 

= j<fxL[g,^. (2.2) 

By using the degrees of freedom of diffeomorphism, we may take the metric tensor as the following 
conformal flat form 


g^.u = (2.3) 

and thus the zweibeins read 

eoo = eii = eoi = eio = 0. (2.4) 

In this article we always use the conformal metric (12.31) as the background space-time. We assume 
that the space-time is flat both for infinite past and future, i.e. 

^Jim^<,„„(x) = ,„„ = ( J “j 

SO that 


lim II(x) = 0. 

t—>-dioo 

Under the background space-time (12.31) the action (12.21) becomes 

5[T,5] = ^ J d^x • 7 V } • (2.5) 

Here e(x) = detei^(T) and thus ^ The action (j2.5l) apparently possesses the 

following Weyl invariance, 

r n{x) n'{x) = n{x) + w(f) 

1 T(f) ^ T'(f) = e-5^(^)T(f) ■ ^ ^ 

Here we assumed the variation of the metric denoted by uj{x) also satisfies 

lim uj(x) = 0. 

t—^zLoo 
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We investigate the Weyl transformation of the Hamiltonian density. The Hamiltonian density is 
given from the action (|2.2I) as 


H{x) 


Too{x) 
do^ 


dL X dL 

-- 


-L 


• T - L 

4:7V Atv 

H[^,g]{x). 


(2.7) 


The form (12.71) transforms under the Weyl transformation (12.61) as 

dvr dvr \ 2 y dvr dvr \ 2 J 

= H. (2.8) 


This shows that H is Weyl invariant. 

We then quantize the system by Ending the solutions for if) and 'll). The equations of motion from 
()2.5j) read 


f do+ idi 0 '\ / e^^ip \ f do + idi 0 \ ( ^ = n 

V 0 5o - )~ ' V 0 9o - iOi y ) “ 

and when imposing the periodic boundary condition 

'p{t,x + 27 r) = 'ip{t,x), 

X + 27r) = Ppit, x), 


the general solution is given by 


ip{x) 

'Pix) 



+ 00 


n=—oo 


g-if2(x) 
^n(t+ix) ’ 



+ 00 

E *7'”' 


n=—oo 


g-|n(x) 


^ +00 
i;\x) = ^ E 

^ n=—oo 
. +00 
^pUx) = E 


(2.9) 


( 2 . 10 ) 


The suffix [TjH] on the expanding coefficients in (|2.in|) which represents the original held T and the 
background metric g^i, = is added for later use. Thus the eigenfunctions, i.e. the basis of the 

Hilbert space read 




1 e 

1 




1 



h^{x) n{t+ix) 
e , 

^fl{x) ^n{t—ix) 


( 2 . 11 ) 


Here the hermitian conjugation f is understood as the complex conjugation in Hilbert space, after 
making Wick rotation r = —it, 

{ 


gm(T±a:)|t _ ^—in{T±x) 
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The inner product in the Hilbert space is defined on the time-slice t as 

{n\m) = J dx (n|x)[Q][n](a;|m) = j> dx ^{x)eQ f!^\x) = 5n,Tn, 
{n\rh) = J dx (h|x)[Q][n](a;|m) = ^ dx ^/gf}[^'^\x)eo f}^\x) = 6n,m- 

Here the zweibein component 


( 2 . 12 ) 


g 0 ^ g-O(x) 


comes from the diagonal component of the matrix 


/ye ° = 


eo 


which makes the inner product covariant. 

We make the canonical second quantization 


x), x' 
x), ip^t, x' 


1 


ye/ 

1 


-.5{x — x'), 


Ve(i 


■.5{x — x') 


and then we have 


= {5L“i,5eS''"»} = 'S» 


others = 0. 


(2.13) 


(2.14) 


We should notice that the commutation relation (j2.14l) does not depend on the conformal factor 
function Q(x) of the metric (j2.3jl . In other word the algebras of the creation and annihilation of the 
particle are all isomorphic, independent of 12(x). 

In the present article we consider the Euclidean space-time, and in particular in the case of the flat 
space-time fi(x) = 0, the energy operator is given by — and the momentum one is —Thus e.g. 
the eigenstate has the energy and momentum eigenvalues E = n and p = =Fn respectively. 

Therefore left mover V’(/ has the dispersion relation E = —p while the right mover V>(x) has E = p. 
From the above consideration, the vacuum state, i.e. Dirac sea is generically defined as 


|sea)[4,,n] = 

n>0 


Q]t 


| 0 ) 

m<0 


Q]t 


| 0 ) 


(2.15) 


such that 


10) = 10) = 0 for all n. 


It should be noticed that the index n in and is the momentum. Needless to say the 

statement “the Dirac sea is the state filled by the negative energy particles” is meaningful only when 
D(x) = 0, i.e. the background space-time is flat. In this case the conhguration of the particles in the 
Dirac sea is depicted by Figure [2l The operators b\^b \^and are the creation- and 

annihilation-operators of the particles as is shown in Table [TJ 

In order to reveal the transformation property of the Dirac sea under the Weyl transformation, we 
consider the functional representation of the algebra (j2.14|) [B]. As is easily understood the commuta¬ 


tion relation (j2.14l) is expressed in terms of the Grassmannian c-numbers and as 

, _ n[^,n] , [M>,0] _ 3 

On — On 5 On — 




On — On 1 


On — 




(2.16) 


(96». 




5 
















left mover 


''V I 


right movers 

n 




p 


p' 


p 


p 


.o' 


p 


o' 


o. 


O. 


■a 


tx 


tx 

IV 


o. 


'Q. 


XX 


P 


Figure 2: The particle distribution of the Dirac sea 


Then the states in (j2.15ji are written as 

empty vacua: |0) ~ 1, |0) ~ 1, 

Dirac sea: |sea)[,jf^Q] ~ n ^2.17) 

n>0, m<0 

In conclusion, the empty vacua |0) and |0) do not include any wave function component, so that it is 
nothing but “classical number” 1. Therefore it does not transform under the Weyl transformation. 

3 Weyl anomaly from Dirac sea 

3.1 Identity for the Weyl transformation 

In this section we investigate the vacuum-to-vacuum amplitude under the Weyl transformation (12.31) 
in order to show that the Weyl anomaly actually stems from pair production from the Dirac sea. For 
later convenience we rescale the complex fermion field T by using the metric tensor in the following: 

The rescaled field T is introduced so that the Dirac sea becomes invariant under diffeomorphisn:13- 
This invariance of the Dirac sea is shown in the Appendix O Then the Weyl invariance (12.6p of the 
action (12.5D leads to 

5[T,5] = ^ j (fx (3.2) 

Originally the rescaled field dSU was introduced by K. Fujikawa in order to make the path integral measure invariant 
under diffeomorphism. 
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spectrum I 

^ create a positive-energy particle, 

annihilate a positive-energy particle. 

spectrum H 

create a positive-energy particle, 

annihilate a positive-energy particle. 

spectrum HI 

^ create a negative-energy particle 

(annihilate a positive-energy hole), 

annihilate a negative-energy particle 

(create a positive-energy hole). 

spectrum IV 

^ create a negative-energy particle 

(annihilate a positive-energy hole), 

annihilate a negative-energy particle 

(create a positive-energy hole), 


Table 1: The interpretation of the creation and annihilation operators 


where 


n(x) —= n(x) + uj{x), 

T(f) ^ T'(f) = (3.3) 

Now, by denoting the time ordering operator by T and using the expression 

(see (12.71) and (12.81) 1. the following identity for the vacuum-to-vacuum amplitude holds, which will be 
shown soon below: 

[§',c'] (sea|Te" ^ |sea) 

= ,o] (sealTe- ^ |sea) 

= f2j(sea| • Te"-/''^""'^['^’3]|sea)[,j (3.4) 

This identity is reminiscent to the formula (7.13) in the Fujikawa and Suzuki’s text [7] for the Fujikawa’s 
method of anomalies. Here we should notice that the second suffixes VL' and ff of |sea)j^, Q,j and 
|sea)[^ o] represent the metrics on which the Fock spaces are constructed. We denote the Fock spaces 
'Hq,! and T-in that are built up from the vacua |sea)j^, and |sea)j,j^j respectively. Needless to say, 
the inner products and the hermitian conjugations are defined in T-L^i and T-Lq, respectively. As will 
be described later in detail, the numerical factor in the third line will be found to be the Weyl 
anomaly itself which emerges, simply saying, from the difference between the Dirac seas |sea)j,|,, 
and |sea)j,jf^j. It is always possible to insert fields in these amplitudes. However it is not essential 
point so that we disregard on this issue. We may easily understand the reason why the identity (13.41) 
holds by the following explanations: 
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(i) The equality between the hrst and second line of (13.4|) 

We may apply with recourse to the isomorphism of the algebras of the creation- and annihilation- 
operators in the sense that the commutation relations (j2.14l) do not depend on the metric. Thus 

we may replace the operators b\^ ^ by In the path integral formalism this operation 

corresponds to the arbitrary change of the naming of the integration variables. In fact the states 
and the operators in each amplitudes can be written as 




' n'l 


n>0 


m<0 


etc., 






n>0 


m<0 




etc.. (I3.4h l 


As for time translation operators, we write them for the Hamiltonian itself (not for the Hamil¬ 
tonian density). 


g-/d2fH[4-',3'] 

where ^ 


where ^ 




H-cx) 




27r 

n=l 
+ ^0 

e ■> , 

1 +°° - - - - - 

T E»(4T'*4T' + + sLT'sLT* 

n=l 


(I33b) 


In this way the Dirac seas and the normal ordered operators are given by the replacement that is 
the consequence of the isomorphism of the algebras (I2.14p . Note that the renormalized vacuum 

energies that are c-numbers fulfill the equality £'q'^ ^ ^ so that they are regularized 

both on the metric g'^y{x) = 

(ii) The equality between the hrst and third line of ()3.4p 

It is proven by using the Weyl transformation (13.31) in the following manner: As was shown in 
(12.8p . the Hamiltonian density is Weyl invariant 

H[^',g']=H[^,g]. 


We may write in the same manner (I3.4bp as an explicit form, 




H-oo 


where = T ^ + 5L*A5L*A') + Ef-°y 


n=l 


dUb) 


Here the metric is changed to g^uix) = according to the Weyl transformation and the 

vacuum energy is the one that is regularized on this metric. On the other hand, the initial 













and final states |sea)j^, ^,j as t —)• ±oo are given by the condition lim^^ioo = 0, and thus 
are trivially equal to that on the flat space-time, 


|sea)j,j,^j^,j — |sea)j^^^j. 

Of course in —oo < t < -|-oo, Dirac sea |sea)j^, needs to be carefully treated. In fact during 
the time development from t = —oo to f = -|-oo the metric may change, so that the hrst line 
may produce some quantity expressed as (Note that this contribution corresponds to the 
Jacobian factor of the path integral measure in Fujikawa’s method). We postpone to give an 
explicit form to the next subsection. There we will obtain by calculating a change of the 
Dirac sea |sea)j^, under the Weyl transformation. 

In the present subsection we make a preparation to the next subsection in the following. We 
divide the amplitude in the first line of ()3.4p in terms of an infinitesimal time interval At. We 
hereby insert the completeness relation 1 = l’^)[§' n'][§' Fock space T-Ln'^ where 

the n-particle state denotes the basis of "Hq/ and the intermediate state on the each 

time-slice. In this case, the initial and final states are vacua = |sea)j,j, so that when 

we perform summation over all the n-particle states on the time-slice at t, the contributing 
intermediate state is only the vacuum one, i.e. q/]( 0| = |sea)j^, Q/](sea|. We may 

explicitly write the vacuum state (= Dirac sea) on a time-slice at the time t as |sea)j^, 

Then the first line of p.4p reads 


= lim 




sea 


t ^ + oo 




(sea|e“ '^['^'’®'^^*|sea) r^, 






-Ao(sea|e ■^['^'’5']^*|sea)r^, 


• • • X 




(sea|e ]^^|sea) 


['I",0'](T-2At) 

[4",C!'](T-HAt) 


X • • • 


[4",0'](i/+At) 


(seaje 




+ CX) 




( 13314 ) 


N=—oo 


It should be noticed that the time axis t is discretized such that 


t = tjsi = NAt, (—oo < N < -|-oo). 

Next the third line of ()3.4D becomes the amplitude on the metric 

ggu{x) = 

due to the Weyl transformation. Therefore we have to change the basis from 


(3.5) 


1 e 


^ (O + u;) 


into 




gn{t±ii) f • The reason of this change of the basis is because we would like to compare the 
creation and annihilation operators of the fields T' and T which act to the same Hilbert space 

f 1 1 

constructed from the operators expanded by the same basis < ^„%±ix) r ■ We then expand the 
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field in the two ways using two different above mentioned basis. 


i,'[D = ^ E ''I*' 


1 


V r> = 


n=—oo 

T +0O 

1 ^ -r>T/'n'l 6 2 


fM(t+ix) Z 

v « = 






= 71^ E ‘I* 

» r) — 


n=—cx) 


(3.6) 


n=—oo 


l(t-ix) ^ z 

^ r) — 


E 


^n{t—ix) ' 


Then subsequently we change the basis of the projection of the intermediate state |sea) (sea 

As for Dirac sea we may write as 


|sea)[^, = Yl ® n 

n>0 m<0 

= n f / |0)<8) JJ ( / dx^'^(T)[Q/](x|m) ) |0) 


n>0 


m<0 


[C 2 '] (^1^) det [Q/] (x H X n ^'t(f)|o)®n^'^(^)io)- 


m<0 


(&) 


Now the projection operator can be written for x = (t, x) and y = {t, y) as 
|sea)[^, j^,][^, f 2 ,](sea| = det [n'](x|n) • det [n'](x|m) • det(A:|y)[n'] • det(r|y)[n'] 

Lx X JLy y ) 

= ® V^l JjTA'(y)(8)(0| • (I33F) 


Note that the expression (|3.4t1) is manifestly independent of the choice of the basis. Thus 


|sea)[^, f^,](sea| = I JJ 6[f'’^'d|0) ® JJ 


n>0 


m<0 


(oi n ® (oi n ss' 


n>0 


m<0 


n ''p'’">'io> ® n “"i“> 11(01 n ® (oi n ‘S' 

n>0 m<0 


O] 


n>0 


m<0 


— lsea)j^, f^j(sea 


da) 


holds (Note the places of the prime). We may interpret |sea)j^, in the right hand side of 
()3.4h[) as the Dirac sea on the curved space-time of the conformal factor Q'(x) in terms of the 
another space-time basis of the conformal factor D(x). In fact we may regard |sea)j^, as the 
pair production state in terms of the “flat” (i.e. the metric fi) space-time. 

Here we summarize the notation of the Dirac seas in Table [2l We recall the Weyl transformation 
of the conformal component of the metric and the fermion field: 


f}(x) —)• Q'(x) = Q(x) + cu(x), 

T(f) -5- ^'(x) = e^“'^^^T(f). (lOl 


The Dirac seas are defined by the various combinations of the two fermion fields and two metrics. 
While we can dehne another Dirac sea |sea)j^ j^,j, we don’t need it in the current article. 
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Dirac sea 

operators and its expansions 

view of field 

n>0 m<0 


^ on gfj,y = 

|sea}[^,_n,] = blP' i+|0} ® ^^|0) 

n>0 rri<0 


T' on qPp = 

n>0 m<0 

+ etc. 

4'' on Qyy = 


Table 2: The notation of the various Dirac seas. 


3.2 Pair production from Dirac sea corresponding to the Weyl transformation 

In the present subsection, we calculate the transformation of the Dirac sea |sea)j^ on the time- 
slice t under the Weyl transformation. Hereafter we assume that the parameter oj{x) in the Weyl 
transformation (13.3p is infinitesimal. We now consider the Dirac sea |sea)j,^, r2'](i) held 'h' in the 

metric and then we make a change of basis such that the state becomes 

|sea)j^, in the metric In this case our concern is how the state |sea)j,j, is 

expressed in terms of the creation and annihilation operators of the field ^ in the metric 
(see eq. (l3.4hD l. 

By making use of the infinitesimal Weyl transformation (13.3p . the amplitude in the third line of 
the identity (13.4|) may be obtained as follows. The fields ijj' and ip' should not be expanded by using 

f 1 1 f 1 1 

the basis < ^n(t±ix) r but rather by using the basis < ^n(t±ix) p because in the third line the 
amplitude is the one in the metric gfj,y by the Weyl transformation. Thus 


ip'ix) 
ip'{x) 


1 

1 


-\-oo 

E 




n=—oo 


g|n(x) 
^n{t—ix) ’ 


(3.7) 


see 


We would like to express the expansion coefficients in (13.7p . and h^n in terms of 

and which become the physical creation and annihilation operators in the flat space-time 


as t ^ 4 - 00 . Here and is obtained via 


ip{x) 

ip{x) 


(3.8) 


The change of the basis and thus the change of the expressions of the creation and annihilation 
operators remind us the Bogolyubov transformation. Then we may expect the appearance of J^j in 
the third line of (j3.4l) . We will show below that there appears pair production states from the Dirac 
sea. 


In the first step we derive the relations of the two expansion coefficients b^n and b^n’^^ 
shoud notice that the coefficients b\^ are the components of the right hand side of p.4hl) . 


. We 
Prom 
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(j3.3p . (j3.7l) and (j3.8l) . we find 




27r 


+CX) 




^Tn(t-\-ix) 


/ ^ ^*n,mUrn i 


dx 


^ ^ 1 + ^u{x) ^ 


2ti 

+00 


+cx> 


,o] 


m=—oo 


gif2(x) 

^m{t—ix) 


= (3.9) 

m=—OQ 

Here we defined a matrix representation of the infinitesimal Weyl transformation which is time de¬ 
pendent as 

dx f 1 


Wn,m{t) ^ Sn,m + j dx i^^ix{x)^ {n\x) [Q][Q]{x\m) = 6n,rn + ^ ^ |icj(x)|e(" 


1 


Wn,m{t) = Sn,m "h / dx <j -Uj{x) {n\x) [Q][Q]{x\m) = 6n,m + f ^ ) 2 ^i^) ( ^ 


dx f 1 


^(n—m)(t—ix) 


(3.10) 


We find that the commutation relations for b^n and b^n are not the ordinary 


ones: 


t k 

-1-00 

E 

k=—cio 


dn,k + ^ I ® 


k=—oo 
{n—k)(t+ix) 


6m,k + ^ \ 


= 6n,m + f^ a;(x)e("-”^)(*+“) 


The similar relation also holds for b^n This u correction to the algebra of the harmonic oscillator 
may cause the Weyl anomaly. In terms of these representation of the Weyl transformation the Dirac 
sea of the right hand side of (I3.4hp reads up to the first order of u 


sea 




n '>p'-""io> ® n 4*'""!») 

n>0 n<0 

/ +0O \ / +00 , \ 

n E io)®n E 


n>0 \m=—ao 


n<0 \m=—oo 


= det 


^n.m 


|E(x)| 


X det 

+ E 


6n,m T 


27r [2 

dx f 1 

^ 12 


- n,m>0 

oo{x) \ e(^-™)d-“) 


n,m>0 




n>0,m>0 


dx f 1 


n>0,m>0 


27r (2 
dx f 1 




n<0 


+ E f 2^ (3,11) 


n>0 
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The reason why the factor in front of |sea)j^ in the third and fourth lines of (j.S.lip becomes 

detlTn,m X det is presented in Appendix C. In this expression we notice that |sea)j^, is a 

superposition as a sum of the three types of states: 


and 


n<0 

n>0 


Here and 1^., n)j^ are defined by 






0 


^1 


^ 7['I',n]t7[<E-,0]t 

— "_i 


^n-1 ^-m ^n+1 


^-n+l^fh ^-n-1 




for n > 0, m > 0, 
for ra > 0, m > 0. 


(3.12) 


Let us investigate in detail of the equation (13.111) . In the terms involving and and 

summing up over m and rh in the second line of ()3.1ip , the contribution coming from the representation 
matrix ()3.inp of the Weyl transformation turns out to become the determinant and as the state it is 
united to the Dirac sea |sea)j^ as the third and fourth lines of (13.111) . As for the terms containing 

b^m<d^ in the second line of (13.111) . they create the pair-produced states, the left mover particles with 
energy m and left mover holes (= anti-particles) with energy n that can be written as 


n<0 


At the same time the terms containing creates the pair-produced states, the right mover particles 

with energy tti and the right mover holes (= anti-particles) with energy h, which may be written as 

n>0 

In summary, the second line of (13.111) is reduced to be a linear combination of the Dirac sea and 
the pair production states expressed as the last four lines of (I3.11j) . Some simple examples of the 
dispersion relations of such pair production states are depicted in Figure [3l From now on, the pair 
production states are represented as 


|pair)[§,o](f) 

in short. Strictly speaking, we can use the word “pair production” only when D = 0, i.e. in the 
flat space-time. Here, we use it symbolically but it expresses the point. By applying the formula 
det ^ of a matrix U to (13.111) , we obtain 


sea) 




exp E / dx ^u{x)(^{n\x)in]{t)lQ]it){x\n) + (n|x)[n](i)[n](i)(x|n)) 


n or n>0 


sea 




(3.13) 
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3)[^,n] ® ns<o 


Figure 3: The examples of the pair production states. The white circles in the lower half planes 
represent the holes. 


Obviously the ket vector ()3.13p belongs to the Fock space Tin and its hermitian conjugation should 
be taken in terms of the inner product of TiQ. The expression (j3.13p is the Dirac sea |sea)j^, at 
the time-slice t arising from the Weyl transformation and the factor in front of |sea)j^ will turn 
out to be a Weyl anomaly. Such a Dirac sea stems from pair production of the particles. Strictly 
speaking the pair production state |pair)j^ does not contribute when taking the expectation value 
at each time interval of the amplitude (|3.4hp in the first order of u and thus the pair production state 
does not directly contribute to the Weyl anomaly. However by considering (13.131) , the pair production 
state I pair) takes some norm from |sea)j^, and create the state |sea)j^ which contributes 

directly to the Weyl anomaly. In this sense we may say that the pair production state is the origin 
of the Weyl anomaly. Furthermore if we take into account of the higher order terms of oj, the state 
including states more than two pairs production may directly contribute to the Weyl anomaly. 

The Weyl anomaly is produced by summing up through all history from the infinite past to infinite 
future as (j3.4dll the exponential factor appearing in front of the Dirac sea |sea)j^ at each time-slice 
t = tN of (13.131) . In (13.131) the summation J2norn>o diverge and its divergence corresponds to 
the ultra-violet one of the space coordinate x. On the other hand the product fijy in (|3.41ip may 
also diverge. As we will see in (|3.17p . this divergence turns out to be that of the integration over 
the “energy” k^, i.e. the ultra-violet divergence of the time coordinate t. We should notice that 
this divergence is natural and plausible because we adopt the general covariant theory. In the next 
subsection we shall indeed perform the derivation of the Weyl anomaly in the above method. 

3.3 Evaluation of and the Weyl anomaly 

In this subsection we compute in (13.41) by the method described in the last part of the previous 
subsection B. Firstly we divide the time direction into an infinitesimal time interval and at each 
time in the amplitude (|3.4tip we apply the changing formula of the basis (|3.4hp of the projections. As is 
described there, is obtained by summing up over all time the contribution (j3.13p and its hermitian 
conjugation coming from the Weyl transformation of the intermediate Dirac sea at each time-slice. 
Then reads 


To; = n expJ2x I dx^u}{tN,x)(^{n\x)in](tj^)[n]itM}ix\n) + 

N=—cio norn>0 J 

(3.14) 
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Here tAr is defined in ()3.5I) . The reader should remember that the product is taken by the step 
At = tjv — tjy-i. As we show explicitly in eq. (I3.17|) . this product with respect to N is exponentiated 
and becomes a summation furthermore take the continuum limit, At —)■ 0, and it becomes 

the integration over the time t, f dt. It should be also noted that the overall factor 2 in the exponent 
comes from the fact that the bra and ket vectors both contribute. We should be aware the fact that 
(13.1411 is non-unitary transformation because of the reality of the exponent which is the consequence 
of the reality of the matrix Wn,m in (13.1011 . This is rather naturally interpreted similar to the case 
of Schwinger mechanism [T] and Hawking radiation, because the non-unitarity generates the pair 
production. On the other hand, as will be shown in Appendix B, the conformal transformation does 
not lead to pair production since the contribution by the transformation is unitary (the exponent is 
pure imaginary). 

In the course of the calculation of we face the divergence problem in the sense of Yl'^=o 1- We 
would like to use the heat kernel regularization to tame this divergence. 

We now consider the differential operator in the action (13.21) 

and solve the eigenvalue problem 

P^PXrix) = Xlxr{x) (3.16) 


with the eigenvalues and a set of complete orthonormal eigenfunctions {xr(^)}- Then the product 
Oat appearing in (13.141) . we exponentiate to make it the summation over N. Thus after all it becomes 
the trace of each time t, 


+ 00 


Juj = exp E E / 


n>0 N=—oo ' 
f +00 . 


n>0 N=—cxi' 


dk^\ x) 

TT V 27r 

At 







(3.17) 


Here 


and 


I denotes the basis of the plane wave with respect to the discretized time t]\f = NAt 
-g continuum limit At —)• 0 of ilP^. We take the continuum limit At —)• 0, and 


the summation YIn time step At join to the integration over t, f dt. Then the exponent of 


“Formally the summation E)iv ''' means the trace ''' kAf}- However according to the Pauli’s argument the 

eigenstate of the time \t) cannot exist in our case of the present article. Therefore we use the plane wave 
instead of |tAr). If we treat the Dirac sea not in 2-dimensional but in 4-dimensional space-time, \t) can exist since the 
energy spectrum is continuous and unbounded. The “normalization” of the basis j- is written as (fco|fco} = 

E)jv • V= E)jv At. When we take the limit At —>■ 0, it becomes {ko\ko) = J dt • e“**’“* = f dt, 

i.e. the volume of the time. 
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(j3.17p becomes the integration over the whole space-time, 


Juj = exp 



d^x Lo{x)e 


ik°t 


x|n)(n|3:)[n]p) + [n]{t){x\n){n\x 


—ikPt 



(3.18) 


Note that it is here nnderstood that the dfe^-integration is cut off at as can be seen from (I3.17P 
explicitly. This means that At is really representing an ultraviolet cut-off. This cut-off interpretation 
of At also explains the mysterious appearance of (13.141) in which the various factors -the exponentials- 
contain no compensating dependence on the stepping distance in time At = tjv — t^r-i, so that 
gets strongly depending on this stepping distance At (Interpreted as cut-off this strange dependence 
may be O.K.). The kets \n) and |h) are the complete set in —oo < n < -|-oo. In fact the quantities 
of (• • •) in the exponent of the first line of (13.181) is independent of n, because [f 7 ](x|n) is just a plane 
wave essentially. Thus we may rewrite 5 J2n=-oo there the completeness relation 

~ Yn ~ changing the basis from < I to {xr{x)}, Ju in (13.181) reads 


Jui = exp 


/ 


+00 


d xuj{x) y xl{x)Xr{x) 


r=0 


(3.19) 


We then next introduce the exponential regulator of A,, that involves the regularization cut-off M. 
The divergent quantity J^j in (I3.19p is now redefined as 


= exp 


/ +00 ^2 

d^x uj{x) xlix)e~^Xr{x) 

r=0 


Finally by expanding in powers of the cut-off and inserting ^ and ^JgR = 

we obtain 


(3.20) 

—2d^d>^Q{x) 


Jui = exp 


= exp 


lim [ d^xoj{x) I- — • 2d^d^^n{x) 

1^00 j y 12 


1 

27r M^oo 

— lim [ d'^x u}(x) < -\ - ^/gR 

2tt m^oc / ^ M 12 ^^ 


(3.21) 


The finite part of (|3.2ip is proportional to the exponentiated Ricci scalar and we may expect that it is 
the Weyl anomaly. However, we need the relation with the trace of the energy-momentum tensor. 
Indeed we relate the Ricci scalar in the exponent of (13.211) and T)) by using the second and third lines 
in the identity (j3.4l) . 

We in fact derive the Weyl anomaly by using the second and third lines of the identity (|3.4I) . 


(sea|r exp ( — j d?x ) |sea) = (sea| • T exp ( — 

Here the limit 


j d^xH[^,g^ |sea). (3.22) 


lim H(T) = lim ooix) = 0 
>±00 £^±00 


for the calculation see e.g. Ref. [7]. 
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is taken. Thus the Dirac sea of the initial and final states in (j.S.22p are those on the flat space-time 


I sea) _ 

The is obtained by renormalizing the divergent term proportional to in (I3.2ip into the cosmo¬ 
logical constant term, 


= exp 



(3.23) 


Remembering the argument in the subsection 13. II the terms in the left and right hands sides are both 
obtained by integrating the same field T. However we rewrite them by using different metrics to 
produce the different effective actions = f dt 


g-5eff[9'=e2(f'+“)] ^ . g-Seff[9=e2«]_ 


(3.24) 


By taking logarithm on both sides, the difference of the effective action S'eff becomes to the differen¬ 
tiation with respect to the infinitesimal cu, 

I (_2a;rp = J . (3.25) 

This is the equality what we are looking for, i.e. the relation of the Weyl anomaly: 

<T) = (3-26) 


4 Conclusion and outlooks 

4.1 Conclusion 

In the present article we presented that when the Dirac sea in the background space-time with the 
metric is expanded in terms of the Fock space basis of the space-time with the metric 

diiu = there appears the state of the sum of the Dirac sea and a state called |pair)j,^ In par¬ 

ticular, in the case of D = 0 , the state |pair)j^ describes the one of the particle-antiparticle pairs. 
In perturbation theory, between the formulation in terms of the empty vacuum and the formulation 
in terms of the Dirac sea in which the negative energy states are all filled, there is no difference in the 
result. However in the case of the nonperturbative phenomena such as the particle-antiparticle pair 
production, there appears apparent difference as we have shown. The origin of this comes from the 
fact that the vacuum is not empty but there is a sea of the infinitely many negative energy particles 
and these are transformed under the Weyl transformation. 

When the Dirac sea, which is defined in a diffeomorphism invariant manner, is acted by the Weyl 
transformation, the pairs of the particle and antiparticle are created and the Weyl symmetry is broken. 
We have shown that this is precisely the origin of the Weyl anomaly. In the Fujikawa’s method based 
on the path integral formalism [7], the path integral measure violates the symmetry and the Weyl 
anomaly stems from the Jacobian of the Weyl transformation. In our operator formalism, the pair 
production from the Dirac sea on each time-slice, which occurs during the transition from —oo-past to 
J-oo-future, breaks the Weyl symmetry. The Weyl anomaly is obtained by summing up the breaking 
effects over all the time. 

In this manner, we show the formulation of the physical relation between the pair production 
state and the appearance of the Weyl anomaly. Our formulation can be applied to many physical 
phenomena, irrelevant to the detail of the models. For instance the state expression of the pair 
production can be used by changing the gravitational field to the electromagnetic field to investigate 
the dynamics of the electrons and holes in the materials. As for the physics of the high energy particle 
and the universe we are under investigation and some of the topics will be shown in the following 
section lOl 
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4.2 Outlooks 


Really the main motivation for the present work is the hope that we may extend our formulation to 
slightly different but related issues: 

(1) The pair production of bosons to generate the Weyl anomaly: 

In our previous works [6], we formulated the boson sea which is the vacuum being filled with 
negative energy bosons. If we set the background metric to the conformally flat one and transform 
the boson sea by the Weyl transformation, we can construct the pair produced state from the 
boson sea vacuum, which would be the bosonic version of (I3.11|) . All we should do may be 
parallel to the method shown in the present article. 

(2) All other anomalies such as the gravitational anomaly being generated via the pair production 
from the sea: 

We can apply our method to the physical interpretations of all anomalies and unify them in 
terms of the pair production from the sea vacuum. 

(3) The Hawking radiation and the gravitational wave: 

When we consider the sea vacuum in the Schwarzschild or de Sitter space-time background, we 
can express the Hawking radiation as the pair production state. In particular, in the case of 
the complex boson and the boson sea, we can represent the gravitational wave by using the pair 
production state as the Hawking radiation, because the number of the degrees of freedom of 
the gravitational wave is two which is identical to that of the complex boson. Then we may be 
able to evaluate the entropies of the black hole and the inflationary universe using the notion of 
the typicality [8] or the formulation of statistical mechanics based on the thermal pure quantum 
states [S]. 
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A Invariance of the Dirac sea under the diffeomorphism 

The rescaling (|3.1I) in Section [3] is meant that Dirac sea becomes invariant under diffeomorphism. In 
this appendix we give a proof of the invariance. 

Under the infinitesimal coordinate transformation 

x'^ = — e^{x) (A.l) 

the rescaled fermion field T = transforms as 

T'(f) = T(T) + (A.2) 
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As is done for (13.91) the transformation rules of the creation and annihilation operators under the 
diffeomorphism (IA.2P is written as 


dx 




+ CX) 


+ 00 


/.fj] 




/ ^ xyn,m^rn J 


dx 




+ 00 


m=—oo 


+ CX) 


^m{t+ix) 


',n] e 2 


^r 2 (x) 


^m(t—ix) 


— -^n,m^rn 

m=—oo 

Here we defined the representation matrices of the infinitesimal diffeomorphism as 

Dn,m = dn,m + dx (n|x)p ^6^(9^ + [a\{x\m), 


(A.3) 


Dn,m — dji^m A ® dx (n,|x)|'Qj ( €^df^ A I [f2](^l^) 


(A.4) 


By using these expressions, we can estimate the transformation law of the Dirac sea under the diffeo¬ 
morphism up to the first order of the infinitesimal parameter e as 


sea 


D 


n>0 n<0 

= n f E io) ® n ( E ) io) 


-1-00 


n>0 \m=—oo 

= det 


n<0 \m=—oo 


dn,m + j)dx [n](f)(a:ln) ("ik)[Q](t) 

dn,fh, + j> dx (h|a:)p(t) + ^df,eA [n](t)(x|m) 


- n,m>0 


X det 




sea 


ri.,m>0 




= exp 


Ylfdx |[f7](t)(a;|n) (^k)[C!]W 


+ (^|a;)[f7](i) [o](t)(a;|n) 




+ |pai")^,Q]W 


(A.5) 


Then, when we consider the vacuum-to-vacuum amplitude, we obtain the following factor for the 
diffeomorphism instead of in (I3.14p for the Weyl transformation. We also change the basis to 
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{Xr(^)} and finally 


JS = exp 


= exp 

= 1 . 


- 1-00 . 

Y^ / (fxdf,^€f^xlix)Xr{x)'^ 

^ _n ^ 


(A.6) 


In this way the contribution from the Dirac sea to the potential anomaly like the identity (13.41) vanishes 
as total derivative. 


B Invariance of the Dirac sea under the conformal transformation 

In this appendix B, we consider the transformation of Dirac sea under the conformal transformation. 
As is explained briefly below, in the conformal field theory, it is a field theory on the flat space-time. 
Thus it may be reasonable to suppose that the conformal transformation does not contribute to the 
pair production unlike to the case of the Weyl transformation in curved background space-time. 

Firstly we make the action (12.21) by diffeomorphism such that the metric tensor becomes the form 
of the conformal flat (12.31) . We then use the complex coordinate {z,z) by making the coordinate 
transformation. 


z = e^ = e^°+“ 


z = e'^ = “ 


(B.l) 


The metric ()2.3[) is assumed to be of the form of the product of holomorphic and anti-holomorphic 
part. 


It is written in terms of the zweibein ej^(x) such that 






= 622 = 0. 


Then the action reads 


S = ^ J dzdz d d 

The action ()B.4p is further reduced to be of the flat space-time one 


— [ i 

i-TT J 


S = — I dzdz 
dvr 




(B.2) 


(B.3) 


(B.4) 


(B.5) 


by applying the Weyl invariance 


U]{z) —)• Oj’{z) = U]{z) -|- ^(z) 
uj{z) — u:'{z) = uj{z) + ^{z) 
ip{z,z) ip'{z,z) = e~^^^^^^'ip{z,z) 

'tjj{z, z) Pp'{z, z) = z) 


(B. 6 ) 


We now construct the Dirac sea by making the standard canonical second quantization. The held 
expansion is performed under the periodic boundary condition 

'ip{t,x + 2 tt) = 'ip{t,x), 
ip(t, X + 27r) = x) 
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and given as 


+CX) 


+CX) 




- ~n+2 

n=—oo ^ 

+ CXD 


n=—OQ 

+CX) 




- ~n +2 

n=—oo ^ 


The commutation relations of the creation and annihilation operators are of the usual one 

l&n, = \bn, Others = 0 . 


(B.7) 


(B.8) 


Thus the Dirac sea is constructed as 


|sea) = 6+10) ® 6l^|0), 

n >0 m <0 

s.t. bn\0) = 6„|0) = 0 for all n. (B.9) 

Now the action (IB.5[) is invariant under the following infinitesimal conformal transformation 


z ^ z' = z + e{z) 
z ^ z' = z + e(z) 

^P{z)^i;\z') = {l-^^e{z)}^|;{z) 

i,iz)^i;'iz') = {l-ldeiz)}i:iz) 


(B.IO) 


This transformation can be rewritten as the expansion coefficients of the fields ■0 and 0. In the same 
manner as the case of the Weyl transformation (13.7|) and (j3.8|) , the fields before and after the conformal 
transformation are Laurent expanded by using the same basis as 

ip^z) = {1 + 9e(z)}“2 0 ( 2 ;) — e{z)d'ijj{z) 

+00 


1 - - e(z)a| 


Z ^+2 


E 




^p'(z) = {1 + de{z)] ijj{z) — e{z)d‘4>{z) 

2 'i 


= <; 1 - ^de{z) - e{z)d \ ^ 


^n+i 


+ CX) 

E 


V. 


- ~n+| ' 

n=—00 ^ 


(B.ll) 
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We then obtain the transformation rule of bn and bn, 


bL 


dz 

27ri 


1 - ^deiz) 




+ 00 


^ ^ Cfl,771^7711 


m=—oo 


+ 00 


+ 00 7 

bm 


m=—oo 


^ ^ Cn^TTibm’ 


m=—oo 


Here the matrix representation of the conformal transformation is defined by 

Cn,m = dn,m " / | ^5e(z) + e(z)a| 

a,™ = dn,m + / + 6-(z)a| 


(B.12) 


(B.13) 


As a byproduct, it is easily shown that (|B.12I) is indeed the Bogolyubov transformation in the following: 




+ O0 




t 

m.k 


k=—oo 


-de{z) + €{z)d 


^—k— 


1 

2 


= 5 


n,m 


= 6 . 


n,m 


= 6 . 


n,m' 


X 



dz 

2iTi 


1 

2 


-de{z) + e{z)d 



“ / ^ I ^9e(z) + e(z)5| 



(B.14) 


Similarly satisfies the same property as (IB.14jl . 

We now consider the transformation of the Dirac sea under conformal transformation. The proce¬ 
dure of the calculation is much the same as that of Weyl transformation. In the vacuum-to-vacuum 
amplitude pair production state does not contribute as is shown in Weyl transformation so that we 
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abbreviate it in the following: 

/ +00 


+CXD 


Isea')*" = 


n>0 \m=—cxD 

= det 


n 1 E i») ® n E i“> 

\m=—oo / n<0 \m=—oo / 

Sn,m - / 

X det 
+ I pair) 


n,m>0 




sea) 


- n,m>0 


= exp 


n>0 




X exp 


+E 

n>0 


2-”-i I ^di(E) + e(2)e 1 2”-i 


I sea) + I pair) 


c 


= exp 


n>0 


dz fl„e( 2 ;) 6(2;) 


Vi — Va 

I 27ri 1 2 2 ; 


— n- 


E i dz f 1 p,e(z) e(z) 

i — < -d^ + n ^ ' 


n>0 


27ri I 2 z 


Z 2 


I sea) + I pair) 


(lBd5k l 


In (IB.lShI) we drop the total derivative term with respect to 2 : and 2 : and then insert the Laurent 
expansion of e and e such that 

e(^) = 

n 

£(f) = 


The equation (jB.15hj) reads 


Isea')*^ = exp 


= exp 


(eo - eo) ^ 


n 


n>0 


I sea) + I pair) 


c 


-Imen 

12 


I sea) + I pair) 


(B.15) 


In the last equality we used the C-function regularization X]„>g ^ (IB.ISP the factor in front 

of the Dirac sea |sea)t is nothing but the phase and we may conclude that the conformal transforma¬ 
tion corresponds to the unitary transformation except for the pair production state. Furthermore in 
the factor eo and eo are the zero modes of de{z) and de{z) respectively. Among them the physical 
meaningful terms are only the sum de{z) + de{z) which is the change of the conformal part of the 
metric. We may set Imeo = 0. Thus under the conformal transformation we conclude that as far as 
we concern the vacuum-to-vacuum amplitude (j3.4p . the Dirac sea does not produce any contribution 
like the pair production. This result coincides with the fact that the conformal field theory is just a 
closed theory in the flat space-time. 


C About the determinants in front of |sea) 


[T,D](t) 


in ( IVTTD 


We start from the second line of (I3.11jl and treat the left and right movers independently. For the left 
mover, we can assume that the creation oeprators of the Dirac sea (see e.g. (I2.15p ) is arranged 
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in the ascending order as to the index n: 


|sea)j^ o](f) = • • • |0) (g) [right mover part]. (C.l) 

Then, taking notice of m > 0 of the summation cutting off the indices n and m up to finite M 
and utilizing the fermionic nature of we get 


M / M 


n 1 E |o) 


n=0 \m=0 


IQ ll 

det 


M \ 

) 






rx,m>0 


(C.2) 


The terms about m < 0 become the states in the fifth line of (13.111) . The same procedure can be 
followed for the righr mover. After taking the limit M —)• +oo, we obtain (13.111) . The determinants in 
(I3.4bp can be obtain in the same way except for the continuum index x. 
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